Ising models and topological codes: classical algorithms and quantum simulation 
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We present an algorithm to approximate partition functions of 3-body classical Ising models on 
two-dimensional lattices of arbitrary genus, in the real-temperature regime. Even though our algo- 
rithm is purely classical, it is designed by exploiting a connection to topological quantum systems, 
namely the color codes. The algorithm performance is exponentially better than other approaches 
which employ mappings between partition functions and quantum state overlaps. In addition, our 
approach gives rise to a protocol for quantum simulation of such Ising models by simply measuring 
local observables on color codes. 



Introduction. — In recent years, several new cross 
connections between classical spin models on lattices and 
quantum information science have been discovered [IHS] • 
This line of research has lead to a transfer of knowl- 
edge and techniques between these two fields, yielding 
e.g. insights into measurement-based quantum computa- 
tion [51 0J [TU] , completeness results for classical models 
[5], but also new classical and quantum algorithms [B}- 

Enm. 

A particularly interesting connection has been estab- 
lished between topological quantum systems and the clas- 
sical Ising model [2-4J. Topological quantum states, 
most prominently Kitaev's toric code 12J, have re- 
ceived tremendous attention recently. They constitute 
new phases of matter [13] and are essential for topo- 
logical quantum computation [14] . The connection be- 
tween topological quantum systems and the classical 
Ising model is obtained by considering the Ising model 
partition function: the latter can be represented as the 
inner product between a topological quantum state and 
a certain product state [2 4J. Here the topological state 
encodes the geometry of the Ising model and the product 
state encodes interaction strengths and temperature. 

In this paper we show how a classical-quantum corre- 
spondence can be exploited to design a new algorithm for 
estimating the partition function of Ising models. More 
precisely we consider Ising models with 3-body inhomo- 
geneous interactions defined on two-dimensional lattices 
embedded in surfaces of arbitrary genus. Such models 
are known to be linked to the topological color codes 
(TCCs) [3]. Even though our partition function algo- 
rithm is purely classical, it is obtained by utilizing the 
associated quantum formulation. The particular connec- 
tion to TCCs exploited in the algorithm is a modification 
of the overlap mapping described above. Interestingly, by 
means of this modification, our algorithm offers an expo- 
nential speedup as compared to algorithms (both quan- 
tum and classical) that are based on using state overlaps 
directly [1 [TTJ ES] (see also [HI HI])- Furthermore, the 
algorithm applies to Ising models in a real-temperature 
regime. This is in contrast to other approaches that only 
work for complex i.e. unphysical temperatures [7J. 

What is more, our approach also leads to an efficient 
quantum simulation of such 3-body Ising models: the 



partition function can be estimated by measuring expec- 
tation values of certain simple local observables of the 
topological color code state. The different parameter 
regimes of the Ising model can be accessed by varying 
the measured local observables, providing an easy tool 
to study phase transitions or other interesting features 
of 3-body Ising models. Interestingly, TCC states can be 
prepared efficiently on a quantum computer (i.e. with 
a polynomial number of gates; however, also any other 
method of preparing such topologically protected systems 
immediately allows for an efficient quantum simulation of 
the corresponding classical spin model). 

We finally remark that both our classical algorithm 
and quantum simulation are obtained by exploiting a par- 
ticular symmetry of the TCC i.e. its self- duality. As such, 
our method is not restricted to color codes, but applies to 
all classical spin models associated with self-dual quan- 
tum codes. Interestingly, the toric code however does not 
fall into this class. 

Color codes. — Here we define TCCs [15) . A 2-colex 
C is a 2-dimensional lattice embedded in a torus of ar- 
bitrary genus g, with the following properties: (a) every 
vertex of the lattice has degree 3, and (b) the faces of 
the lattice are 3-colorable. Examples are the hexagonal 
lattice (Fig. [I]) and the 4-8 lattice. 
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Figure 1: (a) Color code on a hexagonal lattice: qubits are 
placed on vertices, and operators Xj and Zf act on all qubits 
of a face, (b) Corresponding Classical Ising model: spins are 
placed on faces, and 3-body interactions take place between 
spins on faces f a , g a and f) a adjacent to vertex a. 

With every 2-colex we associate a Hilbert space by 
placing a qubit on each vertex. Denoting by X and Z 
the standard Pauli matrices, for every face f we define 



the following commuting operators (FigflJ): 



Xt '.= 



I X v and Zi := 



> Z v . 



(1) 



vef 



vef 



The associated TCC is the space of all states \ip) satis- 



fying X 



ZM) for all faces f. Equivalently, 



this is the ground space of the Hamiltonian 



n TCC = -J2 ( x f + z f)- 



(2) 



f 



It will be important in this work that this Hamiltonian 
is self-dual in the sense that "H T cc is invariant under the 
transformation A — > HAH applied to each qubit sepa- 
rately, where H denotes the Hadamard operation (this 
follows from the property HXH = Z). In general, the 
ground state space is degenerate. In this work we are 
interested in one particular ground state, namely 



|fi) = TV" 



"IK*" 

f 



*f)|0) 



(3) 



where N > is a normalization factor. 

3-Body Ising models. — Consider an arbitrary 2- 
colex with vertex set V and face set F. Place a classical 
Ising spin er^ G {1,-1} on every face f G F . At each ver- 
tex, we consider a 3-body Ising interaction between the 
3 spins located on the faces adjacent to each vertex (Fig. 
[I]). More precisely, for every vertex a G V let f a , g a and 
f) Q denote the 3 faces adjacent to a (labeled in no partic- 
ular order). Then the spins af a , a Ba and <7f, interact as 
— Ja&fa&Ba&i)* where J a is a (positive or negative) site- 
dependent interaction strength. Altogether, the energy 
of a spin configuration a = {erj : f G F} is 



#0) = ~Y1 JaO U o% a o ha . 
aev 



(4) 



The partition function of this model is Z = ^ ct e~^ H ^ a \ 
and all relevant quantities of the classical model can be 
derived from Z. 

The connection between TCC states and the above 
Ising model is obtained as follows. Define 



pPJa. 



\a u 






with a G V 



aev 



(5) 



With the product state \a) :— ® a \a a ) we then have 

Z = <y(Sl\a). (6) 

Thus the partition function is related (up to the easily 
computable prefactor 7) to the overlap between |0) and 
the product state \a), which contains information about 
temperature and couplings. The mapping (pi) was first 
proved in ^j in a slightly different form, i.e. using a differ- 
ent product state instead of \a), in which case it applies 



to surfaces of trivial topology (g = 0) . The identity ffib 
however holds for arbitrary genus [18] (Sec. [Til ). 

The overlap mapping (|6| can be used to compute the 
partition function Z by calculating the overlap (0|a). 
In fact, one can (approximately) compute such overlaps 
either on a quantum computer or on a classical com- 
puter. A natural quantum algorithm approach is to com- 
pute (0|a) via the Hadamard test. Since TCC states 
are stabilizer states on V qubits, there exists a quan- 
tum circuit C of size poly(V) such that \il) = C\0) P^] . 
Since \a) is a product state, there also exists a poly-sized 
quantum circuit C such that \a) — C'\0). This leads 
to (Q\a) = (0|C+C|0). Since the overall circuit OC has 
polynomial size, the matrix element (0|C^C|0) can be es- 
timated on a quantum computer using the Hadamard 
test. This allows to estimate (^|a) in polynomial time 
with error l/poly(y) (with a success probability which is 
exponentially close to 1). Using Z = 7- (0|a), this leads 
to an estimate of Z with error 



iold 



poIy(V) ■ 



(7) 



An equivalent approach was taken in [5] for Ising models 
associated with the toric code state. 

The quantity (0|a) may however also be estimated di- 
rectly on a classical computer. This can be done using 
the probabilistic techniques from [llj . where a general 
method was given to estimate overlaps between stabilizer 
states and product states. It turns out that the resulting 
classical algorithm allows to estimate Z in polynomial 
time with the same error scaling (|7| as achieved by the 
Hadamard test. 

In the present paper we will provide a new classical 
algorithm to estimate Z which outperforms the above 
approach by an exponential factor, as shown next. 

Main results. — Let V denote the number of vertices 
of the colex and F the number of faces. We provide an 
algorithm with runtime 



0(V 3 -\-l og^-) 

£ A 1 — p 



(8) 



which outputs an estimate Z EST of the partition function 
Z, such that the inequality 



Z\< 



7 



VW^ 



(9) 



is satisfied with probability at least p. The algorithm 
applies to all temperature regimes and arbitrary inho- 
mogeneous couplings J a . 

We now discuss the performance of our algorithm and 
compare it to other approaches. First, the 3-body Ising 
models considered in this work are exactly solvable (us- 
ing the Bethe Ansatz) for hexagonal and 4-8 lattices on 
trivial topologies in the case of uniform couplings |20j . 
Our result, however, is considerably more general in that 
it applies to arbitrary lattices on surfaces of arbitrary 
genus and for site-dependent couplings. 



Second, the runtime scaling (fsl) implies that in poly(F) 
time it is possible to achieve e = l/poly(V) and a success 
probability p which is exponentially (in V) close to 1. The 
error in estimating Z as given in pi thus scales as 



A = 



y/2 F ~ 2 • poly(F) ' 



(10) 



We can now compare the performance of our algorithm 
with approaches that estimate state overlaps directly, 
giving rise to the approximation scaling (17b. Interest- 
ingly, the accuracy of our algorithm constitutes an im- 
provement by a factor of 2~ F ' 2+1 compared to the over- 
lap approach. Note that this factor is exponentially small 
in the number of sites V of the lattice in many cases of 
interest. Indeed, it can be shown [TS] (Sec. VI I that 



F= V —^L + 2 



(11) 



where g is the genus of the surface. Thus whenever the 
genus is not too large (e.g. g constant, g = 0(\ogV) or 
even any g < V/5) we will have F = 0(V), in which case 
2-F/2+1 w jjj k e ex p 0nen tially small and the error A is 
thus exponentially smaller than A id- 

What is more, we find [TS] (Sec. VI) that 



Z < 



7 



vw^ 



(12) 



Hence, whenever 2~ F / 2+1 is exponentially small, the ap- 
proximation error A,^ is exponentially larger than Z 
itself and thus meaningless. This demonstrates that, for 
such lattices, the overlap approach is not useful for any 
values of temperature and couplings. In short, we have 
demonstrated that our algorithm yields drastic improve- 
ments over the overlap approach. 

An important question is to understand for which 
regimes of temperature and couplings the partition func- 
tion Z is comparable in size to A. In particular, if A 
turns out to be larger than Z, the approximation pro- 
vided by the algorithm is not useful. Note that it is hard 
to determine whether or not such an issue occurs, since 
computing Z is hard — indeed this is the goal of the al- 
gorithm in the first place. This type of problem is not so 
much a drawback of our algorithm in particular, but is 
typical for all algorithms which, as ours, provide so-called 
additive approximations (see [7]). Here we show that, 
for certain instances, A is indeed provably significantly 
smaller than Z. An example of such a case is T = and 
Ja > 0. This is the zero temperature regime in a fer- 
romagnetic system. In this case one has Z = r y/2 F / 2 ~ 1 
[TB] (Sec. VII) which is indeed much larger than A. A 
e is T = oo. In this case it can be shown 
(Sec. |VrJ) that Z = 1 2- F / 2+1 4r9. This is much 



second exa mp 



larger than A as long as the genus g scales as 0(\ogV). 
In addition to the above examples, further instances of 
meaningful approximations can be identified. In partic- 
ular, configurations of the parameters /3 and { J a } which 
are sufficiently close to the aforementioned examples will 



also give rise to a meaningful approximation as well, us- 
ing a continuity argument. 

Proof ingredients. — The proof of our main result 
is given in [18] (Sees. [HI IV). Here we outline the main 
steps. Step 1. First we will prove a new identity re- 
lating the overlap (fi|a) to a certain quantum expecta- 
tion value. More precisely, using the shorthand notation 
\oc a ) = x a \0) +y a \l) where \a a ) is as in (pi), we introduce 
the tensor product operator 



A:-- 



)A n 



where A„ 






We will show that 



V2 F ~ 2 • (fi|a) = (n\A\n). 



(13) 



(14) 



This identity is proved by exploiting the self-duality of 
color codes. This self-duality implies that the state |fi) 
is an equal superposition state of the form 



|n)oe5>) 



(15) 



ses 



where, crucially, the set S C 1\ is a classical self- 
orthogonal linear code. It is the self-orthogonality of S 
which will allow to relate the overlap (fi|a) to the expec- 
tation value (Q|A|f2). Interestingly, this seems to be a 
specific feature of color codes. 

Combining ( 14 ) and (pi) yields a relation between the 
partition function and the expectation value (fi|vl|fi): 



Z = 



7 



v/2 



F-2 



(n\A\n). 



(16) 



Note that this identity has a different character from the 
overlap mapping ([6]). Indeed, writing out the overlap 
(f2|a) = 5Z(f2|x)(x|a) yields a direct term-by-term cor- 
respondence to the partition function Z in the sense that 
each term (f2|a;)(a;|a) is positive and immediately identi- 
fied a Boltzmann weight of the corresponding Ising model 
(up to the multiplicative factor 7). In contrast, the ex- 
pansion (fi|A|n) = ^(f2|x)(y|J7)(x|A|y) contains both 
positive an negative terms. Thus a nontrivial recombi- 
nation of terms occurs to ensure that the resulting sum 
indeed yields a proper partition function. 

Step 2. We will show that the quantum expectation 
value (fi|j4|fi) can be approximated efficiently on a clas- 
sical computer with error e with runtime (|8| . To achieve 
this, we will use methods for classically simulating quan- 
tum systems. A key ingredient is that each matrix A ai 
which is a real orthogonal matrix with determinant — 1, 
can be decomposed as A a = ZU^D a U\ here Z is the stan- 
dard Pauli matrix, U = HP where H is the Hadamard 
gate and P = diag(l, i), and D a is a diagonal matrix. As 
a result, we find 



(a\A\Q) = {Q,\Z® V [U 



®V\TJ®V^ , 



D a U® v \Q). (17) 



Furthermore Z® v \tt) = |fi) so that the operator Z® v 
can be absorbed. Writing \tp) := U® v \fl) and f(x) := 
(a;|0D o |x) yields 



(Q\A\Q) = J2m\{^)\ 2 = {f)- 



(18) 



This is the expectation value of a random variable on the 
set of bit strings x which takes the value f(x) with prob- 
ability |(x|<£i)| 2 . The expectation value (/) can hence be 
estimated by sampling this random value; the number of 
samples needed to get an error e is 0(l/e 2 ). Apart from 
the number of samples, one needs to take into account 
the complexity of generating each sample as well as the 
complexity of computing f(x). The following property 
is now crucial: since U is a Clifford operation and the 
TCC states are stabilizer states, also \ip) is a stabilizer 
state. The Gottesman-Knill theorem then allows to sam- 
ple the distribution {|(x|<p)| 2 } efficiently, more precisely 
in 0(V 3 ) time [21] . Furthermore, its easy to show that 
the values f(x) can be computed efficiently (more pre- 
cisely: in 0(V) time). Putting together the time com- 
plexity of the sampling, the time complexity of comput- 
ing f(x) and the number of samples yields a classical 
algorithm with scaling (ISj) to provide an e-approximation 
of (fi|j4|fi). Using (16) then also yields a classical algo- 
rithm for estimating Z with error Q. 

Quantum Simulation. — The identity ( 16 ) also gives 



rise to an immediate method for a quantum simulation 
algorithm of the classical 3-body Ising models, i.e. to 
estimate Z. The simulation simply consists of prepar- 
ing |fi) and measuring the expectation value of the local 
observable A which encodes coupling strengths and tem- 
perature of the model. The state |fi) can be generated 
efficiently on a quantum computer. Since TCC states are 
stabilizer states, a poly-sized quantum circuit with 0(V 3 ) 
gates suffices [19] (however, any other method to prepare 
color codes, e.g. as ground states of some effective Hamil- 
tonian, is also suitable). Furthermore A can be measured 
in 0(V) time — in fact all individual observables A a can 
be measured simultaneously. Using standard probability 
theory bounds, this results in a quantum simulation that 
allows one to estimate (Q|A|fi) with probability p and 
error e with runtime (ISj) . Using (16) then immediately 



yields a quantum simulation algorithm for estimating Z 
with runtime (pi). Thus the quantum simulation per- 
formance is the same as that of the classical algorithm. 
Nevertheless, the quantum simulation might be appeal- 
ing in its own right, e.g. owing to the simplicity of the 
protocol. 

Conclusion and outlook. — We have presented a 
classical algorithm for the simulation of 3-body Ising 
models associated with topological quantum systems. 
The algorithm was constructed via a detour, i.e. by clas- 
sically simulating properties of quantum systems. It is 
capable of simulating inhomogeneous models on lattices 
with arbitrary genus, and can be applied to systems in 
the real-temperature regime. This opens the way to effi- 
ciently simulate such systems and investigate their prop- 
erties. At the same time, our approach gives rise to 
an efficient quantum simulation of the classical model 
that only involves measurement of local observables on 
a topological color code state. We have shown that the 
achieved approximation (both for the classical algorithm 
and quantum simulation) is meaningful in the low and 
high temperature limit. Furthermore our method gives 
rise to an exponential improvement as compared to pre- 
vious approaches. It is interesting to note that our tech- 
niques do not apply to 2-body Ising models, which are 
associated with Kitaev's toric code via an analogous over- 
lap mapping. Despite the fact that the overlap mapping 
is very similar, the underlying classical code is in the 
case of the toric code not self-orthogonal which prevents 
the application of the techniques established in this pa- 
per, i.e. rewriting state overlaps as expectation values. 
The same fact also prevents a direct quantum simulation. 
Whether our methods can be generalized to other models 
remains an open problem. 

Finally, we mention that the techniques presented in 
this work (in particular Eq. ( 14)) can be used in a rather 
different area, namely to compute the geometric measure 
of entanglement in TCC states. This is done elsewhere 
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Supplementary Material 

Here we provide mathematical proofs of our results. The 
structure of this Supplementary Material is as follows. In 
section III we introduce preliminary concepts. In section 
ll| we prove the overlap mapping Eq. J6J). In sections 
III| and |IV| we prove our main results (Step 1 and Step 
2, respectively). In section W\ we provide our quantum 
simulation algorithm. In section we prove Eqs. (Ill and 
Eq. (12). Finally, in section |VIIJ we discuss the high and 



low temperature behavior of the 3-body Ising partition 
function as discussed in the last paragraph of the Main 
Results section. 



I. PRELIMINARY CONCEPTS 

In this section we introduce the basic concepts that 
will be central in the proofs of our results. After hxing 
some basic notation in section |I A] we introduce topolog- 

we 



ical color codes (TCCs) in section IB In section IC 
recall basic properties of Calderbank-Shor-Steane (CSS) 
codes, of which TCCs are examples. 



A. Notation 

We will denote by A and Z the standard Pauli ma- 
trices. If u = (ui, . . . ,u n ) £ Z 2 is an n-bit string, we 
will often consider the following n-qubit tensor product 
operators: 



)X U " 

Z u n 



X(u) := X U1 ®- 

Z{u) := Z Ul ®--- Z ,[ " (19) 

where X 1 — X and A = / and similar for Z. We call 
X(u) an A-type operator and Z(u) a Z-type operator, 
respectively. It is easily verified that, for every u, v G Z 2 , 
we have 

X(u)\v) = \v + u) 

z(u)\ V ) = Mr T » 

X(u)Z(v) = {-l) uTv Z(v)X(u) (20) 

where \v) denotes an n-qubit standard basis state in the 
usual sense and where u + v is computed over Z 2 . 

B. Topological color codes 

Recall the definition of a 2-colex C given in the main 
text. Let V and T denote the sets of faces and vertices 
and F and V denote the number of faces and vertices, 
respectively. The following two properties easily follow 
from the definition of 2-colexes: 

Lemma 1. (Basic properties of 2-colexes) Consider 
an arbitrary 2-colex. Then every face contains an even 
number of vertices. Furthermore, every two distinct faces 
overlap in precisely two vertices, or do not overlap at all. 



Consider the TCC associated with the colex C. The 
operators Xf and Zf (recall (fTl) in the main text) are 
called (A-type and Z-type, resp.) face operators. It is 
known that there are precisely F — 2 independent A- 
type face operators (and hence F — 2 independent Z- 
type operators as well) [15] . The dimension of the TCC 
(number of encoded qubits) is 2 m where 

m = # qubits — # independent stabilizers (21) 
= V-2(F-2) (22) 

= 4 - 2 X = Ag, (23) 

where x = F + V — E is the Euler characteristic of the 
surface and where in the last identity we have used that 
X = 2 - 2g. 

In this work we will be interested in the state |f2) as 
given in ([3]) in the main text. This state is one of the 
ground states of % TCC i.e. one has 



x f \n) = |n) = z f \n) 



(24) 



for all faces f. In the following it will also be convenient 
to work with the unnormalized state |fi') := JV|f2). 

Next we derive an explicit expansion of |f2) in the com- 
putational basis. To do so, define the V x F vertex-face 
incidence matrix B, which describes the incidence rela- 
tion between faces and vertices of the colex, as follows: 
the rows of B are indexed by vertices, the columns are 
indexed by the faces; the matrix element B v t is defined 
by 



B 



»>f 



1 if v € f 

otherwise. 



(25) 



Equivalently, each column of B is a 0/1 vector where the 
Is occur to those vertices which comprise a given face. 
Let Zf denote the set of 0/1 vectors with entries labeled 
by the faces: t = (tf : f E F) with tf <G Z 2 . With these 
definitions and using (20), we find 





f tezf f 




Y[X?\0) = \Bt) 
f 


hence 





(26) 



(27) 



|JV> = H(I + X f )\0) 



= EIWioHEi^)- 



Finally, the set 



S :={Bt : teZ$} 



(28) 



(29) 



is a linear subspace of Z^: for every s,s' e S it follows 
that s + s' G S where the sum s + s' is computed over 



1,2 ■ Since there are F — 2 independent face operators 
Xf , there are F — 2 linearly independent columns in the 
matrix B. Hence the cardinality of S is |5| = 2 F ~ 2 . 
Putting everything together, we arrive at the following 
identities; 



N 2 = (Q'\n f ) = 2 F+2 



(30) 
(31) 

(32) 



C. CSS codes and states 

A linear subspace C C Zj is called a (classical) binary 
linear code of length n. The elements of C are called its 
codewords. The orthogonal complement of C is the set 



C x := {»€Z": u T v = for all u £ C} 



(33) 



which is also a linear code. A Calderbank-Shor-Steane 
(CSS) quantum code is defined as follows [H]. Let C 
and D be two binary linear codes of length n such that 
D C C . The CSS quantum code associated to the pair 
(C,D) is the space of all n-qubit states \i/j) satisfying 

X(u)Z(v)\ip) = \ip) for every u £ C, v £ D. (34) 

The set of operators 

S:={X{u)Z(v):u£C,v£D} (35) 

is a commuting group called the stabilizer of the code. A 
generating set of S is obtained as follows. Let 



{u 1 ,. •-,«*} C C and {v\...y} C D 
be bases of C and D, respectively, and denote 
o-. t :=A(u l ) and r,- := Z(v J ). 



(37) 



II. PROOF OF THE OVERLAP MAPPING ^ 

Recall the definition of the vertex- face incidence matrix 
B and consider a vector Bt where t — (if : f G F). The 
vector Bi has one entry per vertex of the colex. For a 
vertex a, let f a , g a and \) a denote the 3 faces adjacent to 
a as described in the main text. Then the a-th entry of 
Bt is given by the sum 

*f. + Ha + Ha ( 4 0) 

Recalling the expression ( [28] ) for |0'), it follows that 

\M) = z2 ( S)l t f« + h a +h«)- ( 41 ) 

t a 

Now define \a' a ) := e^ J "|0) a + e~ /5j£ '|l) a . We claim that 



z = (n'| 



(42) 



where Z is the partition function of the 3-body Ising 
model defined on C as defined in the main text. To prove 
this overlap relation, consider an arbitrary 0/1 vector 
i £ 1,2 with entries indexed by the faces of the colex. 
Such a i corresponds directly to a configuration of the 
Ising spins; in particular, define <7f := (— l) tf i.e. if = 
iff <7f = 1. With these notations, we have 



^U +ha +h*Wa) 



aP J o 



if U + i„ 



-f* J ' if if„ + t B 



exp[j3J a a fa a ga a tla 



+ ha=0 
+ ha =1 

(43) 



(36) This implies that 



]J(ha +ha +haK) =Y[cxp[/3J a a u a Sa ai ta }. (44) 



Then the operators {a%, . . . , Cfc,ri, . . . ,rj} form a gener- The identity |42|) follows, 
ating set of S. 

For every binary linear code C of length n, define an erly normalized states. For every vertex a define 
n-qubit state 



Finally, we rewrite the identity ( 42 1 in terms of prop- 



\C):- 



1 



= Y\u) 



(38) 



IHOII Ve 2 ^ + e- 2 ^ 



(45) 



Any state of this kind is called a CSS state. One can 
show that, for every u £ C and v £ C x , we have 



and denote \a) — (££) \a a ). Recall the identity (30) relat 



ing the unnormalized state |fi') with the properly nor- 
malized state 1 51). Together with (42) we find: 



X(u)Z(v)\C) = \C). 



(39) 



Furthermore \C) is the unique state (up to a global phase) 
satisfying the equations (39). Thus a CSS state is a one- 
dimensional CSS code. 

Finally, coming back to color codes, it follows from the 



Theorem 1. (Overlap mapping) Consider an arbi- 
trary TCC on a 2- colex. Then 



discussion in section IB that every TCC is a CSS code 
associated with the classical codes C = S = D, where 



where 



S was defined in (29). Identity (32) shows that |Q) is a 
CSS state; more precisely, we have|il) = \S). 



z = 7 .(f>H 



a£V 



2/3 J a 



(46) 



(47) 



Theorem [T] is related to a similar result obtained in 
[4j . There the authors obtain a relation of the form Z — 
5 ■ (Cl\a) where 5 is an easy-to-compute prefactor as 7 
in theorem nl and where \a) — (^) \a a ) is a product state 
where 



\a a ) = cosh/3J a |0) + sinh/3J a |l) 



(48) 



Note the difference between the product states |a) and 
I a). In contrast to theorem 111 which holds for arbitrary 
2-colexes, the relation Z = ■ (fl|a) obtained in [3] only 
holds for those colexes for which the state |f2) is the 
unique ground state. This means that the dimension of 
the code is 1, corresponding to g = or, equivalently, 
\ = 2 owing to (23). The relation between Ref. [I] and 
theoremfTlis obtained as follows. If the TCC has |fi) as its 
unique ground state, the self-duality of the TCC Hamil- 
tonian T-i TCC implies that the state |fi) must be self-dual 
as well, i.e. H® v \Vl) =\Q). Noting further that 



H\a a ) oc \a a 



(49) 



the identity ( 46 1 implies that 

Z = 1 -(n\a) = (n\H® v \a)=6-(n\&) (50) 

for some S. 



have the following general result. Consider a CSS code 
with stabilizer generators <j\, . . . , a^, n, . . . , tj where each 
Oi is an X-type operator and each Tj is a Z-type operator. 
Define the code Hamiltonian 



n ( 



and define the state 



J2 ai 



E-i 



iv^>=iv- i na+^)io) 



(51) 



(52) 



in analogy with the TCC state |fi). Then the following 
holds: 

Lemma 3. (Self-duality and self-orthogonality) If 

^code is self-dual then \ip) is a CSS state where the as- 
sociated classical code is self-orthogonal. 



Proof: write <r, = X(u l ) for some u l £ IQ. Since 
the code Hamiltonian is self-dual, it follows that k — 
I and (after possible relabeling of the stabilizers) r, = 
Z(u l ). Let C C Zj be the code generated by the u l . 
Using an argument analogous to the derivation of ( 32 ) 
one shows that \ip) — \C). Since Ti\ip) = \i/j), it follows 
that Z(u l )\C) = \C). This last identity implies (using 
Q) that 



III. PROOF OF MAIN RESULT: STEP 1 



Next we prove identity (16) in the main text. To do 
so, in section III A we first show that TCC states are 



CSS states for which the underlying classical codes are 
self-orthogonal. This property will then be used to prove 
(16) in section IIIB 



A. TCCs and self-orthogonal classical codes 

Here we show that the TCC state |fi) = \S) is a CSS 
state of a special kind: the code S is self-orthogonal. This 
means that every two codewords are orthogonal: s t = 
for every s,t E S. Equivalently, S C S . 



Lemma 2. (Self-orthogonality) The code (29) is self- 
orthogonal. 

Proof: to show that S is self-orthogonal it suffices to 
show that B T B = 0, that is, every two columns of B 
are orthogonal and every column is orthogonal to itself 
over Z2. But this immediately follows from lemma [T] 
first, since every face has an even number of vertices, 
it follows that each column of B is orthogonal to itself; 
second, since every two distinct faces either overlap in 
two sites or do not overlap at all, it follows that every 
two distinct columns of B are orthogonal. □ 

The self-orthogonality of S is closely related to the self- 
duality of the color code Hamiltonian H TCC . Indeed we 



nee 



(-ir u \v) 



v£C 



\v). 



(53) 



It follows that v T u l = for all v 6 C and for all i. Since 
the u 1 generate the code C, it follows that v T u = for 
all u,v G C. Hence C is self-orthogonal. □ 

Finally, we will need the following basic property. 

Lemma 4. Let C be a self- orthogonal binary linear code 
of length n. Then Z® n \C) =\C). 

Proof: let d G Z 2 l denote the all-ones vector. Then 



( 20 1 implies that 



Z® n \u) = (-l) dT » 



(54) 



for every u £ Z£. Note that x 2 = x for every x £ Z 2 . 
Therefore 



d T u 



E M * = E 



2 T 

u 4 = u u. 



(55) 



□ 



If u £ C, the self-orthogonality of C implies that u T u — 
0. Hence d T u = for every u £ C. Together with (54) 
this implies Z® n \C) =\C). 



B. The expectation value mapping 



Next we prove identity (16). Consider an n-qubit sys- 
tem and let \n) = (^) |7Tj) be an n-qubit (real or complex) 



9 



product state, where |7Tj) = a^O) +bj|l). Define an asso- 
ciated n-qubit tensor product operator as follows: 



A:= 



i A 



where A := 



a l hi 
b, -a, 



(56) 



Now consider an arbitrary self-orthogonal linear code C 
of length n and the associated n-qubit CSS state \C). 
The following lemma relates the overlap (C\ir) to the 
expectation value (C|A|C): 

Lemma 5. (Overlap rewriting) Let C be a binary 
self-orthogonal linear code of length n and let \n) be an 
n-qubit product state. Then 



'\C\ ■ (C\n) = (C\A\C). (57) 

Proof: For any u = (ui, ... ,u n ) £l>2 we denote 

a a b u = H a, H bj. (58) 

i:Ui— j:Uj=\ 



Noting that ZAi = a^ + biZX, we have 
Z® n A= J2 a a b u X(u)Z(u). 



uezii 



We now claim that 

(C\X(u)Z(u)\C) 



1 if u e C 
otherwise. 



(59) 



(60) 



To prove the claim, first consider u G C. T hen u € 
C 1 - since C is self-orthogonal. It follows from (391 that 



X(u)Z{u)\C) = \C). Second, consider u £ C. Then for 
every v € C we have u + v £ C as well. It follows that 
the state 

X{u)Z(u)\C) = -J= £(-l)« T > + l ,) (61) 

vl c l „ ec 

is a superposition where each basis state |u + u) corre- 
sponds to a bit string lying outside of C . This proves the 
claim. 

Since C is self-orthogonal, we have Z® n \C) = |C) ow- 
ing to lemma |I] With (J59| and (J60| it follows that 



(C\A\C) = {C\Z® n A\C) 

= ^ a fi 6"(C|X(u)Z( M )|C) 



One the other hand, directly applying the definition of 
\C) we find 






This proves the result. 



(63) 

D 



We now return to the TCC state \D) = \S). Consider 
the product state \a) = \a a ) as in (45 1. The associated 



tensor product operator A = A a is given explicitly by 



A 



v^Wl 



-2,3 J a 



o0Ja 

~HJ a 



-13 J a 
P f3Ja 



(64) 



Recalling that |5| = 2 F 2 (see section IB), theorem [II 
and lemma [5] immediately imply: 

Theorem 2. (Expectation value mapping) Consider 
an arbitrary TCC on a 2-colex. Then 



Z = 



vw^ 



(n\A\n). 



(65) 



IV. PROOF OF MAIN RESULT: STEP 2 

In this section we present our classical algorithm to es- 
timate the partition function of the 3-body Ising model. 
The algorithm is based on the expectation value mapping 
( 65 1 . The algorithm is described in section IV B before- 
hand, in section |IV A| we recall some standard theory of 
Monte Carlo sampling. 



A. Chernoff-Hoeffding bound and Sampling 

The Chernoff-Hoeffding bound is a tool to bound 
how accurately the expectation value of a random vari- 
able may be approximated using sample averages. Let 
X\, . . . Xk be i.i.d. real- valued random variables with 
expectation value E := EJQ and \Xi\ < 1 for every 
i = 1, . . . , K . Then the Chernoff-Hoeffding bound states 
that 



Prob 



1 K 



< e > > 1 - 2e" 



(66) 



We will use the Chernoff-Hoeffding bound in the follow- 
ing setting. Let {p x : x € Zj } be a probability distribu- 
tion on the set of bit strings. Let / be a real function on 
Z2 with |/(x)| < 1 and consider the expectation value 



</>=5>x/0r). 



(67) 



Our goal is to obtain an approximation c of (/) by sam- 
pling. To do so, one follows the following standard pro- 
cedure: first sample K times from {p x }, yielding K bit 



(62) strings x , • • • , x ; then compute the number 



fix 1 ) 



f{x 



K\ 



K 



(68) 



The Chernoff bound guarantees that c is e-close to (/) 
with probability at least p provided that the number of 
samples scales as 



jir = o(4-iog-!-). 

e z 1 — p 



(69) 
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If the computational cost of a single sampling run of {p x } 
is denoted by N samp and if the computational cost of 
computing f(x) on input of x is denoted by AT comp , then 
the total cost of computing c is 



[^samp "T -^ *compJ 



1 



0(^-log 



1 



1 



:)■ 



(70) 



B. The algorithm 

Throughout this section we will consider an arbitrary 
2-colex and its associated Ising model with partition 
function Z. As before, we let V and J- denote the vertex 
and face set, and V and F denote the number of vertices 
and faces of the colex, respectively. The associated Ising 
model is described in terms of the following data, which 
are considered to be given as inputs: 

(a) The V x F vertex-face incidence matrix B; 

(b) the inverse temperature /3; 

(c) the couplings {J a }, which may be site-dependent. 

Additional parameters of the algorithm will be the error 
e of the approximation of Z and the probability p with 
which the algorithm succeeds. 

Next we present our classical algorithm for estimating 
Z. Crucial to our algorithm will be a matrix decomposi- 
tion which relates the matrices A a , which are real orthog- 
onal matrices, to the (single-qubit) Clifford group. Recall 
that the latter is the group generated by the Hadamard 
gate H and the 5 -gate by P :— diag(l, i). We will subse- 
quently exploit that the Clifford group acts in a transpar- 
ent way on TCC states (because the later are stabilizer 
states) to arrive at our classical algorithm to estimate Z. 

Lemma 6. (Decomposing orthogonal matrices) 

Consider a 2x2 real orthogonal matrix O. If det(0) = 1 
there exists a (complex) diagonal unitary matrix D such 
that 



O = P^HDHP. 



(71) 



If det(0) = — 1 there exists a (complex) diagonal unitary 
matrix D such that 



O = ZP^HDHP. 



(72) 



Proof: every O <E 0(2) has one of the two following 
forms, depending on whether its determinant is 1 or —1, 
respectively: 



cos 9 sin 9 
— sin 9 cos 9 



= o 4 



cos 9 sin 9 
sin 9 — cos i 



= O. 



Defining D := diag(e 
straightforwardly. 



), the lemma is verified 

□ 



Here it is interesting to mention that an analogue of 
the above lemma does not hold for general U(2) matrices 
i.e. the fact that O is a real orthogonal matrix is crucial. 



Theorem 3. (Classical algorithm for Z) There exists 
a probabilistic classical algorithm with runtime 



Q(V 3 



1 



log 



1 



1 



;) 



P 



(73) 



which outputs an estimate Z EST of the partition function 
Z , such that the inequality 



-^F,ST Z <* 



7 



(74) 



is satisfied with probability at least p. 

Proof: we use lemma [6J Note from (64 1 that each A a 
is an orthogonal matrix with determinant — 1. Writing 
U := HP, lemma [6] shows that there exists a diagonal 
operation 



D a := diag( e ~ lf \e^) 



(75) 



such that A a = ZW D a U for every vertex a 6 V. Since 
TCC states are CSS states defined by self-orthogonal 
classical codes (recall lemma [2]), lemma El implies that 
Z® v \£l) = |fi). Using this last identity together with 
A a = ZWD a U, and denoting 



it follows that 

(n\A\n) = {n\z® v A\n) = {<p\\ 

Denoting 



(76) 



lA»- (77) 



fi(x) 



Re 



\D a \x) 



have 



f 2 (x) = Im(x\(g)D a \x) 

Px = \(x\v)\ 2 , (78) 

(Q\A\Q) = Y.P*h( x ) +i Y,P*hW 

= </i)+;</ 2 ). (79) 

Since U is a Clifford operation and |fi) is a stabilizer 
state, also \<p) is a stabilizer state. It follows from the 
Gottesman-Knill theorem [21] that the distribution {p x } 
can be sampled in A samp = 0(V 3 ) time on a classical 
computer. Further, its is easily verified that the cost 
of computing fi(x) and fi{x) is N corap = 0(V). The 
sampling scheme described in section |IV A| thus allows 
us to compute numbers C\ and c 2 which are |-close to 
(/i) and ($2), respectively, with probability p, with a 
runtime scaling 



0(V 3 



1 , 1 



P 



(80) 



Then c := Ci + ic 2 is e-close to (fi|A|fi) = (/i) + i(/ 2 ) 
with probability p. Using the expectation value mapping 
( 65 ) it follows that the number 



7 



72^2 



(81) 
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satisfies ( 74 1 with probability p, as desired. 



applying the operation O a can be done in 0(V) time. 
□ Furthermore, the cost -/V comp of computing /(#) (on input 



Note that, interestingly, even though the aim of the 
algorithm is to estimate a partition function i.e. a sum 
of positive contributions, the algorithm takes a "detour" 
via the complex numbers; this is done by introducing the 
diagonal matrices D a . 



V. QUANTUM SIMULATION 

Our quantum simulation algorithm for estimating Z 
will be very simple: it will essentially consist of mea- 
suring the expectation value (fi|.A|fi) by preparing \Q) 
with a quantum circuit and subsequently measuring the 
operators A a . 

Theorem 4. (Quantum algorithm for Z) There ex- 
ists a quantum simulation algorithm to estimate Z with 
the same performance | 73)- (74) as the classical algorithm 
in theorem^ 



Proof: Recall the formula (64 1 for the matrices A a . 
It is easily verified that each A a is a real, symmetric, 
traceless and orthogonal matrix. Hence its eigenvalues 
arc +1 and —1 and there exists a real orthogonal matrix 
O a such that A a = O^ZO a - Computing the V matrices 
O a requires 0(V) time resources. Denoting 



yields 



|0:=®O o |Sl) 



(n|A|n> = m® v \o. 



(82) 



(83) 



The quantum algorithm now basically consists of prepar- 
ing the state |£) and measuring the observable Z® v . 
More precisely, writing 



/(*)=!] ("I)*" and p x = \{x\Cif 



(84) 



aev 



where x — (x a : a G V) denotes a bit string, we have 

(t\Z 9V \Q=J2P*K*) = (f)- ( 85 ) 

X 

The cost -/V samp of sampling the distribution {p x } on a 
quantum computer is 0(V 3 ): indeed, since |0) is a V- 
qubit stabilizer state (as are all CSS states), there exists 
an 0(V 3 ) quantum circuit generating this state. Second, 



of a;) is easily shown to be 0(V). From this point on, 
the proof is finished by straightforwardly applying the 
Chernoff bound, analogous to the proof of theorem [3J □ 



VI. PROOF OF EQS. (|llj) AND <[12J) 

First, the identity 



V - 4g 
F = — — * + 2. (86) 



follows immediately from ( 23 1 



Second, theorem [2] shows that 



Z < 



vw^ 



(87) 



where we have used that |(fi|A|fi)| < \\A\\ < 1 since A is 
an orthogonal matrix. 



VII. HIGH- AND LOW-TEMPERATURE 
BEHAVIOR OF Z 

Here we compute Z in two extremal regimes: zero tem- 
perature in a ferromagnetic system and infinite tempera- 
ture. See the last paragraph in the section Main results. 

Consider first the zero temperature regime in a ferro- 
magnetic system i.e. T = and J a > 0, corresponding 
to A n = Z for all vertices a. In this case one has 



(q\a\q) = (n\z® v \n) = i 



\Q) owing to lemma 



where we have used that Z® v \£l) 
@ It follows that Z = 7/2 F / 2 - 1 . 

The high-temperature regime T = oo corresponds 
to A a — H where H is the Hadamard gate. Recall- 
ing definition ( [29] ) of the classical code S and denoting 
-\T}] one has 



71 l 



(n\A\n) 



(n\H® v \n) = y/\s\-(n\ 

| 5 | 2 -V/2 = 2 F-2-V/2 = 4 _ fl 



®V 



(89) 
(90) 



where: in the second identity we used lemma [5j in the 
third identity we used that -\/]S ; }|f2) — V \s) where the 
sum ranges over all s G S, owing to ( |32[ ); in the fo urth 
2 F - 2 ; in the fifth, we used ((23 1. 



identity we used | yS' | 
In conclusion, (fi|v4.|fi) 



2 ; in the fifth, 
4" 9 so that Z 



